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Particles in turbulence live complicated lives. It is nonetheless sometimes possible to find order 
in this complexity. It was proposed in [Falkovich et ah, Phys. Rev. Lett. 110, 214502 (2013)] that 
pairs of Lagrangian tracers at small scales, in an incompressible isotropic turbulent flow, have a 
statistical conservation law. More specifically, in a d-dimensional flow the distance R{t) between 
two neutrally buoyant particles, raised to the power —d and averaged over velocity realizations, 
remains at all times equal to the initial, fixed, separation raised to the same power. In this work we 
present evidence from direct numerical simulations of two and three dimensional turbulence for this 
conservation. In both cases the conservation is lost when particles exit the linear flow regime. In 
2D we show that, as an extension of the conservation law, a Evans-Cohen-Morriss/Gallavotti-Cohen 
type fluctuation relation exists. We also analyse data from a 3D laboratory experiment [Liberzon 
et ah, Physica D 241, 208 (2012)], finding that although it probes small scales they are not in the 
smooth regime. Thus instead of we look for a similar, power-law-in-separation conservation 

law. We show that the existence of an initially slowly varying function of this form can be predicted 
but that it does not turn into a conservation law. We suggest that the conservation of 
demonstrated here, can be used as a check of isotropy, incompressibility and flow dimensionality in 
numerical and laboratory experiments that focus on small scales. 


Introduction 

The word turbulence is often used as a synonym 
for turmoil and disorder. Inherently, particles in 
a turbulent flow perform an irregular, complicated 
motion. It could therefore come as a surprise that a 
quantity depending on the separation between such 
particles could remain constant during their move¬ 
ment. Of course, due to the chaotic nature of the 
flow, one can expect only a statistical conservation 
of this type, recovered after averaging over veloc¬ 
ity realizations. The subject of the present paper 
is the verification, on the basis of both numerical 
and experimental data of Lagrangian tracers, of the 
conservation law first predicted in [i|. 

Lagrangian conservation laws in turbulence have 
been studied previously and have provided much in¬ 
sight into the breaking of scale invariance in such 
flows (see 0 for a review). However, analytical ex¬ 
pressions for them could be derived only for short 
correlated velocity fields (i.e the Kraichnan model), 
and although they were deduced (numerically) and 
observed in a Navier-Stokes turbulent flow, [3|, they 
were asymptotic laws, holding only when the initial 
separation between particles is forgotten. 

In [l[ (i?(t)“'^) was introduced as an all time 
conservation law for a d dimensional flow, where 
R{t) = |R(t)| is the time t magnitude of the rela¬ 


tive separation between a pair of particles starting 
at a fixed distance Rq. The conservation is expected 
for an isotropic flow and for separations where the 
velocity difference u scales linearly with the distance: 
u (X R, which we will refer to here as a linear flow. 
Previously it was believed that it is only an asymp¬ 
totic in time conservation law and thus would be dif¬ 
ficult to observe [2] (see also 0 for a historic review). 
The result in [l| opens up the possibility to observe it 
and subsequently use it as a check of isotropy and/or 
incompressibility, as well as flow dimensionality, in 
experiments probing very small scales. Physically 
relevant situations for these observations are phe¬ 
nomena occuring around or below the Kolmogorov 
scale, such as, for example, tracer dynamics in cloud 
physics and Lagrangian statistics in the direct cas¬ 
cade of two-dimensional turbulence. 

The invariance of (i?(t)“'^) under the time evolu¬ 
tion can be traced back to a reometrical property of 
an incompressible linear flow [^. For such a flow, for 
each velocity realization, the volume of an infinites¬ 
imal d-dimensional hyper-spherical sector, with ra¬ 
dius i?o and differential solid angle dUlo, at f = 0, is 
gual to its transformation under the flow at time t 

RodQo = R'^{t)dnt. ( 1 ) 

Note that it is due to the linearity of the transforma- 


tion that a spherical sector is transformed to another 
spherical sector. The conservation law directly fol¬ 
lows: for a pair of particles in the flow, starting with 
the separation vector Rq 

(Rit)-^) f {R{t)-^) dno f dnt 

J R^‘^ Sd J Sa ^ ^ 

with Sd the volume of the d — 1 -dimensional unit 
sphere and dflg parametrizing the direction of Ro- 
The first equality is a consequence of the assumption 
that the flow is statistically isotropic, the average 
over velocity realizations for a scalar quantity thus 
being independent of the direction of Ro- In the 
second equality the average and the integration are 
interchanged and equation o is used. 

Recasting the conservation law in the form 
(e“^) = 1 with W = liRR'^/Rg) brings to mind 
the Jarzynski equality [a, |7|, in which W is related 
to entropy production in an out of equilibrium sys¬ 
tem. This hints at the possible presence of a Evans- 
Cohen-Morriss Q or Gallavotti-Cohen Q type fluc¬ 
tuation relation. Below we demonstrate that such 
an extension of the conservation law is indeed possi¬ 
ble for symplectic flows. In particular, since a (lin¬ 
earised) two dimensional incompressible flow is al¬ 
ready symplectic, the existence of a fluctuation rela¬ 
tion is guaranteed and no further assumption, such 
as time reversibility, is required. Resemblance to the 
Jarzynski equality also provides a different perspec¬ 
tive on the latter- it can be thought of as a statistical 
conservation law. 

In the following we provide a direct confirmation 
of conservation of (R{t)~‘^') in an isotropic fully de¬ 
veloped turbulent flow. We first present the re¬ 
sults from a direct numerical simulation of a 2D 
flow, with large scale forcing and friction. For scales 
much smaller than the forcing, in the direct cascade, 
the flow is linear and the theory applies. We show 
that as long as the separation between particles re¬ 
mains in this regime, remains constant as 

well. In addition, we find the symmetry relation 
((i?(t)/i?o)^) = (^{R(t)/ Rq)^^ with g = —(7 — 2 as a 

generalization of the conservation of At 

long times this yields a Gallavotti-Cohen type fluc¬ 
tuation relation as discussed above. 

Next we consider data from a direct numerical 
simulation of 3-dimensional turbulence and study 
particles with initial separations in the dissipative 
range, where approximately oc R^. We observe 
the conservation of ending in a regime 

where no longer converges. At those times 


the separations distribution at small scales is in¬ 
duced by the dynamics in the inertial range rather 
than the dissipative range of scales. 

Finally, we turn to data of SD turbulence from a 
laboratory experiment. Initial particle separations 
are taken at scales comparable with the Kolmogorv 
scale of the flow. We hnd that, in distinction from 
the numerical simulation, the scaling of the velocity 
is far from the approximation (u^) oc R^ for such 
separations. Also, deviations from an isotropic flow 
are observed. Thus, we check instead if an analogue 
of in the form of a power law, exists for the 

scales available in the experiment. We show that a 
slowly varying function of this form can be predicted 
using the statistics of the pairs relative velocity and 
acceleration at t = 0. However, it does not appear 
to be a true conservation law. 


Conservation in 2D -DNS results 

The motion of Lagrangian tracers is numerically 
integrated in the direct cascade of two-dimensional 
turbulence with friction, in a square box with side 
length L = 27r. The flow is generated by a 
large scale, J-correlated, random forcing at a wave- 
number corresponding to the box scale which injects 
enstrophy at a rate Oj. The parameters of the sim¬ 
ulation are taken from 0 - The linear friction is 
sufficiently strong to generate a velocity field with 
energy spectrum exponent close to 4.5 and a power- 
law decaying enstrophy flux, i.e no logarithmic cor¬ 
rections to the leading order scaling are to be ex¬ 
pected 0 , 0 - The viscous enstrophy dissipation 
rate 0 ;^, together with the kinematic viscosity v, de- 
hne the smallest dissipative scale 77 ~ 
which is used as a reference scale. Time is made 
dimensionless with the vorticity characteristic time 
{Z represents the mean enstrophy) 
and consequently the reference velocity is Urj = ij/Trj. 

Fig. [T] shows the longitudinal structure function 
of the velocity, where we denote by u; — u ■ R the 
longitudinal velocity difference at scale R. A scaling 
consistent with a linear flow behaviour is observed 
for R < lOry. 

In stationary conditions, Np = 16384 particle 
pairs are introduced into the flow with homogeneous 
distribution and initial separation Rq and their tra¬ 
jectories are evolved in time. The moments of sep¬ 
aration (i?(f)'3) are computed by averaging over Np 
and Nrun = 100 independent runs. 
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FIG. 1. Second order longitudinal structure function of 
velocity in the 2D direct cascade simulation. The line 
represents the behaviour (uf(i?)) oc 
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FIG. 2. Time evolution of different moments of relative 
separation R{t) as a function of time in the 2D sim¬ 
ulations (i?o = 1-lr;). Different point styles represent 
q = —1.0, q — —1.8, q = —2.0 and q = —2.2 (from 
bottom to top). 

In Fig. [2] we present the time evolution of 
for q = —2.2,—2,—1.8,—1. While the moment for 
q = —2.2 {q = —1.8) grows (decays), doing so expo¬ 
nentially for t > 2Tri, the moment with q = —2 = —d 
is conserved up to a time t ~ lOr,,. At this time the 
average particle separation reaches i? ^ lOiy where 
in Fig. [I] deviates from a linear flow. We ob¬ 
serve that the exponential growth of (R) is inhibited 
at approximately the same time (not shown). On 
the other hand, as can be seen in Fig. [2l the ex¬ 
ponential decrease of which is less sensitive 

to pairs with large separations, lasts throughout the 
observation time. 



FIG. 3. L{q) obtained from fitting {R{t)'^) with 
in the time interval 2 < t/r,, < 10 for the 2D simulations 
{Ro = 1.1??). In order to show the symmetry with respect 
to q = —1, both L{q) (red circles) and L(—2 — q) (blue 
squares) are plotted. 

The exponential time dependence of {R‘^{t)) ob¬ 
served in Fig. [2] is expected to be a long time 
feature of any linear flow with temporal correla¬ 
tions decaying fast enough i, [H. Specifically, 
E{q,t) = In ((i?(t)/i?o)'^), the cumulant generating 
function of /3 = \n{R{t)/Ro), should take at long 
times the form E{q,t) = L{q)t. We remark that the 
finite-time Lyapunov exponent a = Pjt tends to the 
Lagrangian Lyapunov exponent. A, in the long time 
limit |14| . 

We obtain L{q) shown in Fig. |3] by fitting 
(i?(t)«) = for times 0.4 < Xt < 2. Evi¬ 

dently L(q) is perfectly symmetric with respect to 
q = —d/2 = —1, i.e L{q) = L{—q — 2), implying in 
particular that L{—2) = 0 as expected. This sym¬ 
metry can be extended to the symmetry E[q, t) = 
E{—q — 2,t) holding at any time t, as we demon¬ 
strate in the bottom figure of Fig. [5] for t — 0.2r^. 
This is a general property of a two-dimensional lin¬ 
ear incompressible and isotropic flow. In fact, this 
follows from the relation 



which, as we will show, holds for every velocity re¬ 
alization. Indeed, using the definition of E{q) one 
obtains E{q,t) = E{—q — 2,t) for an isotropic flow 
by taking the average of (l3|) over velocity realiza¬ 
tions, conditioned that the averages exist. 

To demonstrate ([3]) we recall that for a linear flow 
with a given velocity realization the following de- 
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composition can be used 


Rjty 


= O' {t)A{t)0{t)Ro 


( 4 ) 


where 0{t) is an orthogonal matrix and A a diagonal 
matrix with entries and 15| !]• For an 


incompressible flow detA = 1 so that in 2D pi = 
—p 2 = p. Using dU the integration over i?o can be 
written as 




d^Q = 


[e^{e)A{t)e{9)Y^^ d0 (5) 


where e = (cos(0), sin(0)) and a change of integra¬ 
tion variable absorbed the additional rotation 0{t) 
in Q. On the other hand, (HD implies 


Ro ) 


fir^n — 




^Rit) 


-q-2 


d^t ( 6 ) 


and Ro/R{t) can be decomposed, similarly to (HJ, 
using the flow backward in time 


R{ty 


= Ri o'i-t)A{-t)Oi-t)Rt. (7) 


The inversion of the linear transformation between 
R(t) and R(0) implies A{—t) = A~^{t) [l^. Then, 
A~^{t) is related to A{t) by conjugation with a rota¬ 
tion matrix since A = diag{e^,e~P). This property 
is actually a consequence of the symplectic structure 
of a 2D incompressible flow- for a symplectic flow 
the eigenvalues of A come in pairs of the form 
0. □ Therefore for any m, shifting the integration 
variable similarly to ([S]) , 


\R{t)) 


dflt = 


[e'"{0)A{t)e{9)] dO 



FIG. 4. Verification of the symmetry relation ®, = 

In {R{t)/Ro), in the 2D simulations. The solid straight 
line corresponds to the analytical prediction. 


The symmetry of E(q) implies a symmetry of 
the probability density function (pdf) P{f3{t)), j3 = 
In (JR{t)/Ro), which can be obtained by first extend¬ 
ing the symmetry of E{q) to complex q and then 
retrieving the pdf by an inverse Laplace transform 
of exp[i?(g)]. Using the symmetry and a change of 
variables in the integration, one gets for any time t 


Rim = y) 

PiPit) = -y) 


= e'y = J{y) 


(9) 


where J(ln (i?(t)/i?o)) is the Jacobian of the trans¬ 
formation © from time zero to time t, in 2D. 

In the asymptotically long-time limit Q, which is 
numerically verified in Fig. |4l turns into an Evans- 
Cohen-Morriss @ or Gallavotti-Cohen Q type fluc¬ 
tuation relation. Indeed, the function L{q) is the 
Legendre transform of the long-time large deviation 
function G(cr) so that the symmetry in L(q) implies: 


which by comparison to m leads to 


G{-a)=G{a)+2a (10) 



Finally, combining equation (jS]) with m = —q — 2 
together with equation (jS]) yields ©■ 


which can also be deduced directly from ®. 
This relation resembles the Evans-Cohen-Morriss- 
Gallavotti-Cohen relation, and even more so the re¬ 
lation found in 1^ , but is different from them. Note 
that we do not assume a time reversible velocity en¬ 
semble and work with an incompressible flow (im- 


^ As this is the basic symmetry leading to the relations m 
and unii, these relations, with the replacement of 2 by 2N, 
would also hold for a linear random/chaotic N-dimensional 
Hamiltonian flow that is statistically isotropic, R denoting 
the separation between two points in phase space. 


^ This can be done since the considerations above did not 
depend on q being real, only on convergence of the integrals. 
It seems reasonable that if they converge for real q they 
would also converge for its complex values. 
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plying zero entropy production in the context of dy¬ 
namical systems). The two dimensionality of the 
flow plays a crucial role in the derivation. 

In Fig. El we show E{q,t) for short times, t < 8r^ 
as a function of q, comparing two distinct initial sep¬ 
arations. For Rq = 1.1 77 all curves cross zero at 
q = —2, up to time t = 8r^, while for Rq = 50?7, 
which is at the border of the linear scaling range, 
the initial crossing is shifted above — 2 , and the in¬ 
tersection point changes with time (see middle Fig¬ 
ure). The appearance of the initial crossing above 
—2 in the latter can be explained using a Taylor ex¬ 
pansion of {{R{t)/ Rq^) around t = 0, see equations 
(IIID,(Iini) and (fHl) in the analysis of the experimental 
results, as well as [ij. 

Conservation in 3D - DNS results 

Lagrangian tracers are introduced into a numeri¬ 
cal simulation of three-dimensional turbulence in a 
cubic box at Re\ = 107. The tracers are placed in 
the flow when it has reached a steady state and their 
trajectories are numerically integrated. Turbulence 
is generated by a large scale, ^-correlated random 
forcing. Small scales are well resolved in the simu¬ 
lation (for which kmaxV = 1-3). More than 130000 
pairs were integrated for times up to about 10 for 
300 realizations each, while more than 30000 pairs 
were integrated for 150 realizations for longer times. 
The two datasets are overlapping for shorter times 
so that consistency of the statistics could be checked. 

In Fig. E] we show the second order longitudi¬ 
nal velocity structure function. Up to separations of 
about 7ri the scaling exponent of (uf{R)) is close to 
2 , although deviations are also observed earlier. 

In Fig. [3 In {{R{t) / Ro^) as a function of q is pre¬ 
sented for three different initial separations and for 
various times. A clear crossing of zero at g = — 3 
can be observed up to time t = 8r,, for Rq = 1.7r]. 
For Rq = 6 . 777 , which lies at the transition from 
the dissipative range, the initial crossing is shifted 
above —3 (see intermediate Figure) a trend which 
is much more pronounced for Rq — 70r] at the iner¬ 
tial range. This is the same trend that appeared in 
the 2D simulation and can be explained in a similar 
way. Evidently, this shifted crossing point cannot 
also correspond to a conservation law, as it is time 
dependent, for both separations. 

Although theoretically (^{R(t) /Rq)~^^ is an all 
time conservation law, in reality after some time 
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FIG. 5. \n {{R{t)/RoY) as a function of q for two dif¬ 
ferent initial separations Ro = I.I 77 (upper plot) and 
Ro = 5 O 77 (middle plot). Different curves represent 
different times. Lower plot: E{q,t) (red circles) and 
E{—2 — q,t) (blue squares) for t = 0.2t,, and Ro — l.lrj. 
Data from 2D simulations. 


pairs reach separations where a linear flow approxi¬ 
mation no longer works, see also Fig. [3 and its dis¬ 
cussion for the 2D flow. We study the exit from the 
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FIG. 6. Second order longitudinal velocity structure 
function for the 3D numerical simulation. The two lines 
represent the smooth and rough behaviors. 


conservation regime of ^(i?(t)/i?o) for i?o = l.Try 
in Fig. m After t ^ 8 — lOr^ wild fluctuations, as 
well as a decrease in value, appear for {{R{t)/RoY) 
with q < —2. This behaviour can be explained by 
the formation of a power-law dependence on sepa¬ 
ration in the separations pdf at small R, see Fig. 
m together with a shift of the average separation to 
larger values. Such a dependence has been shown 
to develop for initial separations both in the iner¬ 
tial range and in the dissipative one in 17|. In Fig. 
IH^a) P{R/Rq) is shown for different times for the 
3D simulations. Indeed the left tail shows a power 
law dependence i?“, with a decreasing until reach¬ 
ing a ^ 2 at times 10 ^ t/r ^ 20. This explains 
why the moment {R~^) starts diverging somewhere 
in this interval. For larger times we observe a further 
decrease in a (not shown), compatible with the val¬ 
ues observed in [3. This power law behaviour may 
be due to the contribution of trajectories which left 
the viscous subrange and returned to R < Rq- In¬ 
deed, P{R/Ro) is expected (at long enough times) 
to have a log-normal shape near its maximum, as 
long as the separation remains in the viscous range 
[^. In our 3D simulations, P{R/Rq) never displays 
a log-normal shape, probably as a consequence of 
the limited extension of the viscous range. On the 
contrary, the PDF’s of separations are clearly log¬ 
normal in the 2D case Fig. |9l(b). Although the 
statistics are not sufficient to distinguish the devel¬ 
opment of a power-law dependence in the left tail in 
the 2D data, the large fluctuations in Fig. [5] point 
towards a behaviour analogous to the 3D case. 
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FIG. 7. \a {{R{t)/RoY) as a function of q for different 
initial separations Ro = 1.7rj (upper plot), Ro = 6.7677 
(middle plot), Ro = 70ri (lower plot). The different 
curves represent different times. Data from 3D simu¬ 
lations. 


Conservation in 3D - experimental results 


In the previous sections we have presented sup¬ 
porting evidence from numerical simulations, both 
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FIG. 8. {{R{t)/RoY) as a function of time for various 
q, the initial separation is Ro = l.Yq. Different symbols 
represent q = —2.6, q = —2.8, q = —3.0 and q = —3.2 
(from bottom to top). Data from 3D simulations. 


in two and three dimensions, for the conservation of 
for an isotropic flow at small enough scales. 
It is interesting to check if similar results can be 
obtained from laboratory experimental data where 
the access to very small initial separation is limited 
by both physical and statistical constraints. It is to 
this end that we turn to data from the experimen¬ 
tal set up described in [18| , where neutrally buoyant 
particles are tracked inside a water tank of dimen¬ 
sions 32 X 32 X 50c?7i^. The turbulent flow, with 
Rex = 84 for our data set, is generated by eight pro¬ 
pellers at the tanks corners and particles are tracked 
in space and time using four CCD cameras. The 
cameras are focused on a small volume of Icm^ to 
resolve the smallest scales (the Kolmogorov scale is 
rj = OAmm). At such small scales, the flow is ex¬ 
pected to be isotropic. We study pairs of particles 
with initial separations of a few rj, for which numer¬ 
ical and previous experimental studies of turbulent 
flows found an approximately linear flow [Bii. 


As isotropy and a linear dependence of velocity 
differences on separation are the two assumptions 
entering the prediction of the conservation of 
the experimental set up described above appears to 
be appropriate to test it. Unfortunately, a direct 
check of these assumptions does not seem to support 
their applicability. 

In Fig. [Tni we present (u^), using the relative 
velocity of particles with separation Rq. Here and 
in the following () denotes an average over pairs in 



R(t)/Ro 


FIG. 9. P{R{t)/Ro) for the smallest available initial 
separations in the 3D (a) and the 2D simulations (b) 
respectively. While the 3D case shows the development 
of power-law tails for small separations, in the 2D case 
all PDFs are close to log-normal curves, marked by solid 
lines in (b) 


an ensemble with a given initial separatiorH. When 
it is applied to functions of the relative velocity or 
acceleration at a given scale it denotes the average 
over the pairs with the corresponding initial sepa¬ 
ration, taken at t = 0. As can be seen from the 
inset of Fig. [TOl if any scaling regime, (u^) oc Rq^ , 
is to be assigned, it would be with scaling exponent 
C 2 = 1-5 rather than (2 = 2. It seems that in or¬ 
der to access the linear flow regime in this system 
even smaller distances should be probed, which is 


® We thank Eldad Afik for suggesting his method of pair 
selection for the ensembles 12^ . which we use in this work. 
We further filter out pairs of particles with a lifetime smaller 
than O.bTri to prevent contamination of the time evolution 
by a change of the particle ensemble 
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FIG. 10. Second order velocity structure function in the 
3D experiment as deduced from the initial relative ve¬ 
locity of pairs of particles with initial separation Ro, 
normalized with the squared Kolmogorov velocity u^. 
Inset: The local slope ciln(u^) /cilni?o as a function of 
the separation. The dashed line represents the scaling 
(u^) oc 


not experimentally accessible at the moment. 

For separations at the transition between inertial 
and dissipative ranges, as those apparently probed 
in this experiment, there is currently no theoretical 
prediction regarding the existence of a conservation 
law. Still, in the spirit of the conservation law in 
the dissipative range, we can look for q such that 
{{R{t)/R q^) = 1. Notice that as a function of q, 
{{R{t)/R q)'^) is convex and at each time t it crosses 
1 at (7 = 0 and either has one more crossing point, 
denoted by q*{t), or none. It would be possible to 
establish the existence of a conservation law, with 
q — q*, a q* exists and is time independent. To get 
an idea of what value q* can take we use the Taylor 
expansion of around t = 0 (no assumption is 

made about the nature of the flow): 


with 


qe{t) = c 


2 (ui) Rq 

t{uf) ’ 


(u^) -h {ai)Ro 

{uD 

( 12 ) 


where ui = u ■ Rq, ai = el - Rq and a is the relative 
acceleration. This implies that the conservation law 
described above can exist only if qe(t) is time inde¬ 


FIG. 11. Average longitudinal relative velocity as a func¬ 
tion of initial separation for pairs of particles at t = 0. 
The velocity is normalized by the initial rms longitudi¬ 
nal relative velocity. Inset: The LHS and the RHS of 
equation (d are displayed as a function of Ro- 


pendent, so that q* = q^ = c. This is of course only 
a necessary but not a sufficient condition. 

The time dependence of qe comes from the first 
order in the Taylor expansion oc {ui). It, as well as 
(a/), should be zero for a statistically homogeneous 
or a statistically isotropic flow. Obviously in an ex¬ 
periment {ui) can never be exactly zero, but it can be 
small enough so that at the shortest times measured 
the time dependent term in q,, would be negligible 
compared to the time independent part, c in (1121) . 
In the inset of Fig. [TT]we present (u/) normalised by 
the rms relative velocity. The strong bias towards 
positive {ui) for Rq < 2ri is probably the result of 
the filtration of pairs with life times smaller than 
0.5r,^. Indeed, pairs with approaching particles at 
small separations are harder to track and are there¬ 
fore frequently short lived as well as more prone to 
errors. A similar problem might also be the cause of 
the positive bias for Rq ^ Sry, where particles with 
a large relative velocity are more common. Then, 
particles approaching each other with a large speed 
reach small scales quickly and are lost more easily. 
When the restriction on the life time of the pairs is 
lifted a negative bias emerges. It therefore appears 
that using (ui) for pairs of particles to determine the 
isotropy or homogeneity of the flow is problematic. 
For our purposes, note that for an observation of the 
conservation of {R’^{t)) to be possible the value of 
(ui) must be as small as possible. Initial separations 
Rq = 2.677 — 377 seem most appropriate. 

For a d-dimensional incompressible flow that is 























statistically isotropic a simpler formula for 
be written. For such a flow (see also 26|) 

which can be used to write 

(u^) dln(u?) 

= id +GiRo))', C2(-Ro) = —7j—^ 
{uf} dinRo 

implying that in (|T^ 

qe = c = 2 — d — C2 {Ro) = —1 — C2{Ro)- 



As a side note, we remark that if C 2 is constant, like 
in the inertial range, the initially slowest changing 
function, R{t)^ = is independent of Ro 

[ij. More generally, requiring only {ui) = 0, the Ro 
independent slowest varying (among twice differen¬ 
tiable functions of R) function F{R), is determined 
by demanding {(PF{R)/dt^') = 0, resulting in 


FIG. 12. The exponent q*{t), deduced by matching 
R{tY ^ = Rq , plotted as a function of Ro at dif¬ 
ferent times. The dashed line represent the prediction 
q*{t) = c, with c from (1121) . Inset: The negative of the 
contribution to c due to {ai), Ro {ai) / (uf), plotted as a 
function of Ro- 


F{R) oc / exp 


-(^t) 

ziuf) 


dz 


dx. 


(16) 


with a an arbitrary separation and = u^—uf. For 
a 3-dimensional isotropic flow this formula reduces 
with the help of (fT51) to 


F{R) oc 


1 


(uf) 


dx. 


(17) 


Here, however, we focus on the simpler functions 
F{R) = which, through qe, are Rq dependent. 

The isotropy assumption is checked using equa¬ 
tion da, as suggested in [26|, in the inset of Fig. 
El where we present both the LHS of this equation 
and its RHS as a function of initial separation. The 
difference between the two is about 1.5 for most sep¬ 
arations, implying that isotropy is violated. It could 
also be that the flow is, in fact, isotropic: although in 
theory measuring relative velocities for pairs of par¬ 
ticles at a given distance is the same as using two 
fixed probes in the flow, in practice it might not be 
so (for example due to correlations entering through 
the method of pair selection). 

In any case, we conclude that the prediction (IT5|) 
should not work for our data, as indeed it does not, 
and instead use in the following the more general def¬ 
inition of c, In Fig. [T2]we present the exponent 
q*{t), as deduced from the data, for five consecutive 
times and compare it to c. As expected, q*{t) ap¬ 
pears approximately constant and equal to c for the 



FIG. 13. The evolution in time of {{R{t)/RoY) for pairs 
with initial separation of Ro — 2 . 677 . The dashed dotted 
curve represents q — c= —1.32 


separations where {ui) is smallest. Note that c has a 
non-negligible contribution coming from the devia¬ 
tion of {ai) from zero. This contribution is displayed 
in the inset of Fig. [121 

We are now in a position to ask if for the sep¬ 
aration where initially q*(t) w c, {R‘^{t)) is just a 
slowly varying function or a true conservation law. 
In Fig. [T31 for the separation Rq — 2.6?7, we compare 
{{R{t)/R oY) as a function of time for different gs, in¬ 
cluding q = c. After a time shorter than {R^{t)) 
noticeably decreases. During this time the average 
separation between particles in a pair changes by less 
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than one percent. It is thus improbable that {R‘^{t)) 
is truly constant up to that time, but rather that it 
changes slowly, meaning it is not a conserved quan¬ 
tity. We conclude that, at the scales accessible to 
this experiment, a power-law-type of conservation 
law probably does not exist but a slowly varying 
function of the separation can be identified. 
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